Magnetic response of nanoscale left-handed metamaterials 
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Using detailed simulations we investigate the magnetic response of metamaterials consisting of 
pairs of parallel slabs or combinations of slabs with wires (including the fishnet design) as the 
length-scale of the structures is reduced from mm to nm. We observe the expected saturation of 
the magnetic resonance frequency when the structure length-scale goes to the sub-micron regime, 
as well as weakening of the effective permeability resonance and reduction of the spectral width of 
the negative permeability region. All these results are explained by using an equivalent resistor- 
inductor-capacitor (RLC) circuit model, taking into account the current-connected kinetic energy of 
the electrons inside the metallic parts through an equivalent inductance, added to the magnetic field 
inductance in the unit-cell. Using this model we derive simple optimization rules for achieving optical 
negative permeability metamaterials of improved performance. Finally, we analyze the magnetic 
response of the fishnet design and we explain its superior performance regarding the high attainable 
magnetic resonance frequency, as well as its inferior performance regarding the width of the negative 
permeability region. 

PACS numbers: 41.20.Jb, 42.70.Qs, 81.05.Xj, 78.67.Pt 



I. INTRODUCTION 

Left-handed metamaterials (LHMs), i.e. artificial com- 
posite structures with overlapping negative permittivity 
and permeability frequency bands giving rise to negative 
index of refraction^, have attracted recently an expo- 
nentially increasing attention. The main reason behind 
this attention is mainly the novel physical phenomena 
associated with those materials (negative refraction, op- 
posite phase and energy velocity, reversed Doppler effect 
etc) , which result to new capabilities in the manipulation 
of electromagnetic waves. Such an important capability 
is the superlensing capability of LHMs^, i.e. the abil- 
ity to offer subwavelength resolution imaging, which can 
have important implications in many scientific, techno- 
logical and every-day life areas, like imaging, microscopy, 
lithography, ultra-compact data storage, etc. 

Since the demonstration of the first left-handed 
material 4 -, in 2000, operating in the microwave regime, 
many left-handed (LH) structures have been created^—, 
and important efforts for the better understanding and 
the optimization of those structures have taken place. 
Among the various efforts within the LHM research, a 
large part has been devoted to the extension of the fre- 
quency of operation of LHMs from the microwaves to the 
optical regime, where the superlensing-bascd applications 
can find an important ground for manifestation. These 
efforts led to metamaterials with negative permeability 
operating in the few THz regime already in 200 4 12 ' 13 , 
which soon were followed by the first structures of neg- 
ative permeability and/or negative index of refraction in 
the telecommunications regime and more recently in the 
lower visible regime. (For reviews of the existing research 



efforts on infrared (IR) and optical metamaterials see 
Refs3!£.) 

While the first and most of the existing microwave 
LHMs are systems made of split-ring resonators^ (SRRs, 
i.e. interrupted metallic rings, giving rise to resonant 
loop-like currents, and thus to resonant permeability 
involving negative permeability values) and continuous 
wires (leading to the negative permittivity response^), 
in most of todays high frequency LHMs the SRRs have 
been replaced by pairs of slabs (or stripes, or wires 
- see Fig. 1(a). Like the SRR, the slab-pair also be- 
haves as a resonant magnetic moment element, where 
the magnetic moment is created by resonant currents, 
antiparallcl in the two slabs of the pair, forming a loop- 
like current. In most of the experimentally realized opti- 
cal slab-pair&wire structures the slabs are as wide as the 
corresponding unit cell side and are physically connected 
with the wires, leading to a design known as fishnet^— 
(see Fig. 1(d)). Fishnet design was able to give the high- 
est in frequency LHMs up to no w 26 ' 27 . 

The main reason behind the replacement of SRR by the 
slab-pair for the high frequency metamaterials, apart of 
slab-pair simplicity in fabrication (which is also a crucial 
parameter), is its ability to exhibit negative permeability 
response for incidence normal to the plain of the pair; 
this makes possible the demonstration of the negative 
permeability response with just a monolayer of slab-pairs. 
Indeed, up to now most of the demonstrated "magnetic" 
metamaterials (i.e. materials of resonant and negative 
permeability) and LHMs arc single layers, while only few 
multilayer samples have been fabricate d 13 ' 28 . (Note that 
in the optical regime what is difficult to be achieved is 
the negative permeability component of a LHM, since 
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the negative permittivity response can be easily obtained 
using metals; that is why most of the existing efforts to 
go to the optical LHMs start from attempts to achieve 
structures of only negative permeability.) 

Since many of the existing attempts to create high fre- 
quency magnetic mctamatcrials and LHMs are based on 
the scaling down of known microwave designs, there are 
already efforts trying to determine the possibilities and 
the limitations of this scaling approach 3 ^ - — . (Note that 
the properties of the metals, which are involved in most 
of todays metamaterials, are drastically different in the 
optical regime compared to microwaves - there, metals 
behave almost as perfect conductors.) Most of those at- 
tempts concern SRR systems and they have led to two 
important conclusion o 30 i 34 ' 35 : 

(a) By scaling down a SRR, the frequency of its reso- 
nant magnetic response does not continuously increase, 
but after some length scale it saturates to a constant 
value. This value was found to be dependent on the 
SRR geometry employed, and with proper modifications 
of this geometry (e.g., adding gaps in the SRR) it could 
go up to a small fraction (e.g. 20%) of the plasma fre- 
quency, i.e. to the middle visible range. The satura- 
tion response of the magnetic resonance frequency was 
explained taking into account the contribution of the ki- 
netic energy of the electrons associated with the current 
inside the SRR ring to the magnetic energy created by 
the loop current (or, equivalently, taking into account the 
dispersive response of the metal in the conductivity). 

(b) The magnetic permeability resonance becomes 
weaker and weaker by going to smaller length-scale SRR 
systems, and below some length scale it ceases to reach 
negative values. This weakening of the permeability was 
attributed to the kinetic energy of the electrons (giving 
rise to saturation) in combination with the increased re- 
sistive losses in the metal as one goes to higher frequen- 
cies; these losses are strengthened by the resonant re- 
sponse, implying long-time interaction of the wave with 
the metallic structures. 

Although the existing studies are very revealing con- 
cerning the high frequency response of magnetic meta- 
materials, they examine the influence of the kinetic en- 
ergy (or the dispersive response of the metals) only to 
the magnetic resonance frequency and not to other fea- 
tures of the resonant magnetic response (like resonance 
shape and damping factor), neither clarified the role of 
the losses in the saturation of the magnetic resonance 
frequency. Moreover, the role of the various geometrical 
parameters in the high frequency response of metama- 
terials still remains to be determined, as to identify the 
dominant parameters determining this response and to 
define optimization rules for those materials. 

In this paper we attempt to study all the above men- 
tioned issues. We will be restricted to systems based on 
pairs of slabs, alone or in combination with continuous 
wires; this is mainly due to the fact that slab-pair-based 
systems arc offered for an easy experimental demonstra- 
tion of negative permeability or negative index response 




FIG. 1: The unit cell of the four designs studied, (a): Narrow- 
slab-pair system; (b) wide-slab-pair system; (c) slabs&wires 
system; (d) fishnet design. Panel (e) is a magnification of 
panel (a) where the structure parameters appearing in the 
structure simulations are shown. The parameters of the pair 
are given as a function of the scale parameter a — (lat- 
tice constant along propagation direction; equal to the system 
thickness): lattice constants oe = 2.97a^, an = 2.19cik, slab- 
length I — 2.19ak, slab- width w = 0.47ak, thickness of the 
metal t m = 0.25cik and thickness of the substrate t = 0.5ak- 
For the fishnet design (d) the width of the slabs is equal to 
the corresponding unit cell side (oh), while the width of the 
"necks" (continuous metallic parts joining the slabs along E 
direction) is w n = 0.469ak- The dielectric spacer separating 
the metallic pairs has been considered as glass (with relative 
permittivity et = 2.14), while for the metal plasma frequency 
and damping factor the aluminum parameters have been em- 
ployed. 



and, moreover, have been proven up to now the most 
promising systems for the achievement of high frequency 
negative permeability and negative index metamaterials. 
The structures discussed here are shown in Fig. 1. Fol- 
lowing the approach of Refsi 30 ' 35 , we will attempt to 
analyze the high frequency magnetic response of those 
structures, to compare their performance and to propose 
optimization rules for them. For that we examine in de- 
tail the scaling behavior of the magnetic resonance fre- 
quency and the magnetic permeability as the structures 
are scaled down from mm to nm scale. 

The basic idea that we will use to reproduce and under- 
stand the small length scale (high frequency) behavior of 
our structures is the consideration of the kinetic energy 
of the current carried electrons. This kinetic energy, as 
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being proportional to the square of the velocity (and thus 
of the frequency, just like the magnetic energy) , is added 
to the magnetic energy of the structures and in small 
length-scales it dominates the magnetic metamaterials 
response. The consideration of this kinetic energy is done 
here through an equivalent "kinetic" inductanc e) 30 ' 36 (or 
electrons' inductance), added to the magnetic held induc- 
tance in an effective resistor-inductor-capacitor (RLC) 
description of the artificial magnetic structures. 

Specifically, the paper is organized as follows: In Sec- 
tion II we present the "high frequency" response of our 
structures, as revealed by numerical simulations concern- 
ing the magnetic resonance frequency, the form of the 
magnetic permeability resonance and the losses. The 
wave propagation characteristics in those structures are 
analyzed and explained in Section III, using an effective 
RLC description of the structures and taking into ac- 
count the dispersive behavior of the metals through the 
kinetic metal inductance. Based on the results and anal- 
ysis of Sections II and III, in Section IV we present basic 
optimization rules for the achievement of high frequency 
magnetic metamaterials and left-handed materials with 
improved performance. There we discuss also the fishnet 
design, which has been proven up to now the optimum 
design for achievement of optical negative index materi- 
als. Finally, in Section V we show that the simple RLC 
circuit model does not only have qualitative power, but 
it can be used also to give quantitative results if plugged 
with accurate relations for the capacitance and the in- 
ductance of the system. 



II. NUMERICAL SIMULATIONS 

In this section we present calculation results concern- 
ing the scaling of the magnetic resonance frequency and 
the magnetic permeability of the structures shown in Fig. 
1. The geometrical parameters used in the simulations 
are those mentioned in Fig. 1; for the permittivity of the 
metal the Drude dispersion model has been employed, 
i.e. e = eo[e^ — Wp/(cj 2 + iwy m )}, with the parameters 

22.43 x 10 15 
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" x , collision frequency j m = 12.18 x 
= 1. Using Drude dispersion model one takes auto 



matically into account the mass and any kinetic energy 
of the current-carrying electrons inside the metal (the 
contribution of the bound electrons is also taken into ac- 
count, through the constant e' )). 

The calculations presented here have been per- 
formed using the Finite Integration Technique, employed 
through the Micro Wave Studio (MWS) commercial soft- 
ware. Using MWS, the transmission and reflection from 
a monolayer of the structure have been obtained; these 
data have been used for the determination of the effective 
permittivity and permeability of the structures, through 
a standard retrieval procedure based on a homogeneous 
effective medium approac h 37 ' 38 . 

In Fig. 2 we present the magnetic resonance frequency 
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FIG. 2: Scaling of the magnetic resonance frequency with the 
linear size of the unit cell along propagation direction (a*,) for 
the four designs of Fig. 1: Fishnet (black circles), wide-slabs 
(red squares), narrow-slabs (green diamonds), slabs&wires 
(blue triangles). The solid symbols indicate the existence of 
negative permittivity values, while the open symbols indicate 
that the permeability resonance is weak and unable to reach 
negative values for the Re(/i). 
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FIG. 3: Frequency dependence of the real (top panel) and 
the imaginary (bottom panel) part of the magnetic perme- 
ability of the wide slab-pair system for various length scales. 
The legends indicate the lattice constant along propagation 
direction for each specific system. 



for the four structures of Fig. 1 as the structures are 
scaled down from mm to nm scale; in Fig. 3 we show the 
real and imaginary part of the magnetic permeability as 
a function of frequency, for wide-slab-pairs (structure of 
Fig. 1(b)) of various length scales. 

As can be seen in Fig. 2, in slab-pair systems we 
observe the same behavior as the one reported earlier 
for SRR structure a 30 ' 35 : while in larger length-scales the 
magnetic resonance frequency scales inversely propor- 
tional to the structure linear size, at frequencies in the 
near infrared (IR) towards optical regime this linear scal- 
ing breaks down, and the magnetic resonance frequency 
saturates to a constant value. This saturation value is dif- 
ferent for the different designs employed, with larger the 
one of the fishnet structure. Note that in the slabs&wires 
case (structure of Fig. 1(c)) the presence of wires does 
not affect the magnetic resonance frequency of the slabs, 
while in the fishnet structure the presence of wires leads 
to higher saturation value for the magnetic resonance fre- 
quency. This behavior of the fishnet design will be dis- 
cussed and explained in Section IV. 

It is important to mention here that the saturation 
values for the magnetic resonance frequency of the slab- 
pair-systems are in all cases larger than the saturation 
values obtained for SRRs of a single gap22 (like, e.g., U- 
shaped SRRs; single-gap SRRs are the only SRR-based 
system that has been fabricated in the nm scale), indi- 
cating once more the suitability of the slab-pair-based 
systems for the achievement of optical magnetic meta- 
matcrials. 

Concerning the permeability results shown in Fig. 3, 
we observe that, just like the SRR s 30 ' 35 , as the length 
scale of the structure becomes smaller the permeability 
resonance becomes weaker, ceasing ultimately to reach 
negative values. This weakening is revealed in both the 
real and imaginary part of the permeability resonance, 
and it will be analyzed in more detail and quantified in 
the following paragraphs. (Note that the "truncation" 
of the resonances of Fig. 3, which is observed in the 
larger scale structures, is a result of the larger influence 
of the periodicity in these length scales^ -. In smaller 
scales this influence becomes smaller, due to the deeper 
sub- wavelength scale of the corresponding structures and 
the weaker resonant response, which results to smaller ef- 
fective index and thus larger wavelength inside the struc- 
tures.) 

One quantity which is of great interest in left-handed 
materials and it is strongly affected by the weakening of 
the permeability resonance is the width of the negative 
permeability regime, which roughly corresponds with the 
operational band width of a LHM. In Fig. 4 we present 
the relative band-width (i.e. band-width divided by the 
lower frequency of the negative permeability band) for 
the four structures shown in Fig. 1. As can be seen 
in Fig. 4, the operational band-width, which in larger 
scales is almost independent of the length scale, in smaller 
scales it is strongly reduced, ultimately going to zero for 
all designs. This shows that the negative permeability 
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FIG. 4: Band-width, Aw, of the negative magnetic perme- 
ability regime vs the inverse unit cell size, l/ak, for fish- 
net (black circles), wide-slab-pairs (red squares), narrow-slab- 
pairs (green diamonds) and slabs&wires (blue triangles). The 
band-width is normalized by the minimum frequency (lower 
limit) of the negative permeability regime, a; m i n . au. is the 
lattice constant along propagation direction. 

is ultimately killed in the nm scale structures. Among 
our four structures the wide-slab-pair one is characterized 
by the larger band-width. Finally, it is worth-noticing 
the reduced band-width of the fishnet design, compared 
to the band-width of the slab-pair-only structure. This 
bandwidth behavior will be discussed in the next section. 

Since a main issue for the achievement of high fre- 
quency metamaterials of satisfactory performance is the 
losses, due to the increased resistive losses in the metal- 
lic components going to the optical regime, we performed 
a detailed analysis of the losses in high frequency meta- 
materials, trying to estimate which aspects of the high 
frequency metamaterial response are mainly affected by 
the resistive losses, and, consequently, to seek ways to 
minimize the influence of those losses, using proper de- 
sign modifications. 

As a first step we calculate the losses as a function of 
frequency for various length scales of our systems. In Fig. 
5 we present these losses for the wide-slab-pair system of 
Fig. 1(b). The losses, A, have been calculated through 
the relation A = 1 — R — T, where R and T are the reflec- 
tion and transmission coefficient, respectively, through 
one unit cell of the structure along propagation direc- 
tion. As is expected, the losses show a dramatic increase 
by going to smaller length scales and higher frequencies. 
This increase seems to have an exponential dependence 
on the magnetic resonance frequency. 

From Figs 2-5, one can see clearly the decreased per- 
formance of our structures going to nm scale, inhibiting 
their ability to give high quality optical left-handed meta- 
materials. An interesting question arising here is the role 
of the resistive losses on this decreased performance. To 
examine this role, we repeated the above shown calcu- 
lations using for the metal a 2tt x 1000 times reduced 
collision frequency compared to that of aluminum, i.e. 
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FIG. 5: Loss, A, per unit cell (%) for the wide-slab-pair sys- 
tem for various length scales of the system, close to the mag- 
netic resonance frequency saturation regime. The legends in- 
dicate the lattice constant along propagation direction. 



7, n = 12.18 x 10 10 /27r. The results concerning the wide- 
slab-pair system are shown in Fig. 6. 

Specifically, Fig. 6(a) shows the saturation of the mag- 
netic resonance frequency and Fig. 6(b) the real part of 
the resonant permeability response for the smaller length 
scales. Comparing the result of Fig. 6 with those for the 
non- reduced value of j m , one can see that the satura- 
tion of the magnetic resonance frequency seems to be 
totally unaffected by the value of the metal collision fre- 
quency. This indicates that the loss-factor of the metal 
employed is not able to affect the highest achievable mag- 
netic resonance frequency of each specific design. On the 
other hand, the metal loss-factor can affect the minimum 
length scale able to give negative permeability response, 
as shown comparing Fig. 6(b) with Fig. 3(a). Indeed, 
in the small "f m cases the magnetic permeability reso- 
nance is quite stronger, maintaining negative values up 
to smaller length scales. Although the strength of the 
resonance (as measured, e.g., by the minimum value of 
the Rc(/i)) seems to be strongly affected by the -f m value, 
calculating the width of the negative permeability regime 
(in the cases that such a regime exists), Aw - not shown 
here -, it is observed that it is only slightly affected by 
the jm value; it tends to zero for both high and low val- 
ues of 7 m , indicating that even in the absence of resistive 
losses one can not go to arbitrarily high-frequency nega- 
tive permeability response over a practical band-width. 

The results presented above raise many questions con- 
cerning the high-frequency magnetic mctamaterial re- 
sponse and the main phenomena and factors determin- 
ing this response. Before attempting an interpretation 
of this response we will summarize here the main effects 
observed so far and the main questions that one needs 
to address as to clarify and to be able to predict the ex- 
istence or performance of optical negative magnetic per- 
meability response: 

(a) The magnetic resonance frequency of slab-pair- 
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FIG. 6: (a): Saturation of the magnetic resonance frequency 
for the wide-slab-pair structure using for the metal a 2000-7T 
times reduced collision frequency compared with that of alu- 
minum, i.e. j m = 7a(/20007t (solid-black line with squares). 
For comparison the corresponding result for 7 m = ja.1 is given 
(dashed-red line with circles). Dotted-dashed line (blue color) 
shows the scaling of the magnetic resonance frequency as it is 
obtained from our RLC circuit model (for the 7ai case), (b): 
Magnetic permeability resonance (shown in Re(/x)) as a func- 
tion of frequency for the wide-slab-pair structure of Fig. 1(b), 
for 7 m = 7^;/20 007T and for three different structure length- 
scales. Legends denote the lattice constant along propagation 
direction, ak. 



based systems, while in mm scale structures scales in- 
versely proportional to the structure length-scale, in sub- 
/jm scale structures saturates to a constant value, de- 
pending mainly on the geometry of the structure and 
independent of the resistive losses in the component ma- 
terials. Among the structures that have been studied 
here the fishnet design leads to the higher saturation fre- 
quency. Moreover, while in the slabs& wires design the 
magnetic resonance frequency is almost unaffected by the 
presence of the wires, it is not the same in the fishnet de- 
sign, where the saturation value is quite higher than the 
saturation value for the the slabs only case. 

(b) The magnetic permeability resonance in sub-^jm 
scale structures becomes more and more weak by reduc- 
ing the structure length-scale, and ultimately ceases to 
reach negative values. The length scale at which /i < 
stops to exist depends on the structure geometry (design) 
and on the resistive losses in the component materials, 
especially in the metallic parts. Among the structures 
that we have examined, negative permeability in smaller 
length scales is achieved in the slab-pair structures, not 
associated though with higher magnetic resonance fre- 
quency. Negative permeability at higher frequencies is 
achieved in the fishnet design. 

(c) Concerning the relative spectral width of the nega- 
tive permeability regime, this width, while it is constant 
in mm scale structures, as one goes towards nm scales 
it becomes smaller and smaller, up to almost vanishing. 
This width depends on the structure geometry and seems 
only slightly dependent on the resistive metallic losses. 
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III. ANALYTICAL MODEL AND 
INTERPRETATION OF THE RESULTS 

To understand and explain the results presented in the 
previous section, we use the common approach of describ- 
ing the artificial magnetic structures (at the resonant 
magnetic response regime) as equivalent effective RLC 
circuits. Using circuit theory and basic electromagnetic 
considerations^, one can easily obtain an expression for 
the frequency dependence of the effective magnetic per- 
meability, fj,(uj), for a lattice of "magnetic" elements^. 

For the explanation of the high frequency magnetic re- 
sponse we follow the approach of Ref^S, based on the 
consideration of the kinetic energy of the electrons in the 
metal, Ek, besides the magnetic energy of the resulting 
electromagnetic field. Thus, we replace the magnetic in- 
ductance of the system in the effective RLC circuit de- 
scription by the total inductance resulting as a sum of the 
magnetic field inductance, L, and the kinetic inductance, 
L e , where L e is defined by E k = N e m e v 2 /2 = L c J 2 /2 (N e 
is the total number of electrons, m e is the electron mass 
and v e the average electron velocity). 

Assuming a slab-pair system like the one shown in Fig. 
1(e), of length I and slabs separation t, excited by a 
magnetic field of the form H = 7Joe _4W *and direction as 
shown in Fig. 1(e), and applying the Kirchhoff voltage 
rule, one can obtain 



(L + L e )I+-I- 



IR = —0 = u) fioltHoe 



(1) 



4> is the external magnetic flux, = fj,oltH, and R and C 
the resistance (frequency independent, accounting for the 
ohmic losses) and the capacitance of the system, respec- 
tively. The obvious solution of Eq. flTJ is I = Ioe~ tut , 
with 



Io = 



uj 2 [nolt/(L + L e )} 



lu 2 -1/(L + L e )C + iuR/(L + L e ) 



(2) 



Having the current one can easily obtain the pair 
magnetic dipolc moment, m = area x current = UI, 
and the magnetization M = (N LC /V)ltI = (l/V uc )UI, 
Nlc is the number of "RLC" circuits in the volume V, 
and V uc = aEflH«k is the volume per unit cell, where 
«E,flH,ttk are the system lattice constants along the E, 
H, and k direction respectively (an > w, oe > I, flk > t) 
- see Fig. 1). 

Finally, using M = Xm(u)H, n{w)/no = 1 + Xm (<*>), 
with Xm the magnetic susceptibility, one obtains that 



, x M (l/Kc)(M0W 2 /(^ + ^)V 2 1 



with 



1 R 



y/{L + L e )& L + L e ' 



(3) 



(4) 



lulc is the magnetic resonance frequency of the system 
and 7 is the dumping factor, representing all the losses 
and the scattering mechanisms. 

The inductance L e can be easily calculated by calculat- 
ing the kinetic energy of the electrons, Ek = N e m e v 2 /2 = 
V w n e m e v 2 /2, and expressing the velocity through the 
current, / = ewt m n e v e [n e is the number density of free 
electrons, e is the electron charge, and V w = wt m l is the 
volume of the metallic slab - see Fig. 1(e)). This way, 
one can obtain L e = lm e /wt m e 2 n e = (l/wt m )(l/u)p£o), 

where uj p = y // e 2 n e /m e Eo is the plasma frequency of the 
bulk metal. 

Note that exactly the same results as in Eq. ([3]) can 
be obtained by considering, instead of the kinetic energy 
for the derivation of L e , the dispersive behavior of the 
metal conductivity. Indeed, starting with the frequency 
dependent Drude-typc conductivity 



ie - 



«7n 



we obtain for the total resistance 



R 



tot 



1 I 

aS 



(- 



s ui 



1 2) c 



R — iu)L e . 



(5) 



(6) 



Using R to t = R — iuL e in Eq. (JTJ in the place of R 
and only the magnetic field inductance (in order to avoid 
counting twice the kinetic inductance L e ), one can obtain 
the same current solution as in Eq. (2). 

In the following, to simplify our discussion we will con- 
sider the magnetic field inductance to be given by the 
inductance of a solenoid of area It and length w (see Fig. 
1), i.e. L = (jftlt/w, and the capacitance by that of a 
parallel plate capacitor of area wl/2, plate separation t 
and dielectric core of relative dielectric constant £&, i.e. 
C = EoEb(wl)/t. (Note that these formulas are appropri- 
ate for the case of wide slabs while describe much less 
satisfactory the narrow slabs case; in any case there is 
a numerical correction factor of the order of one - see 
also Section V. Note also that in our discussion wc will 
not consider any inter-unit-cell capacitance^!, which is 
important in the case of quite long slabs, i.e. slabs that 
approach the unit cell boundaries along E direction.) 

With the above considerations, the equation <j3j) for the 
magnetic permeability takes the form 



with 
F' 



F- 



F'lo 2 



L + L, 



Itw inter-pair volume 
V uc unit cell volume 



(7) 



(8) 



Using the above formulas and their behavior going to 
small length scales we will show in the following that one 
can reproduce and explain all the high frequency mag- 
netic response of artificial magnetic structures. 
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For that, it is important to notice that by scaling-down 
the structures uniformly, i.e. all the lengths scale propor- 
tionally to a basic length a (= here), both the capac- 
itance and the magnetic inductance scale proportionally 
to a, while the kinetic inductance and the resistance scale 
proportionally to 1/a, i.e. 



C — £o£b- 



wl 



a, L 



a 



fJ-o- 



R 



— , and L e 
a 



I 



1 



1 



wt m ui.Eo a 



(9) 



The above formulas shows the increasingly pronounced 
role that kinetic inductance (also the resistance) plays in 
the smaller scales. Specifically, one can see that the ratio 
L/L e is of the order of 40ii m /Ap, with the typical value 
of X p (A p = 2t:c/ijJ p ) being around 85 nm for Al and 130 



for Agii. Thus, for 



smaller than 100 nm the kinetic 



inductance, L e , becomes appreciable and may dominate 
as the length scale becomes smaller and smaller. 



A. Magnetic resonance frequency 

Taking into account the expression for the magnetic 
resonance frequency in Eq. in combination with 

Eqs. ©, it can be shown that the magnetic resonance 
frequency of a slab-pair system has the following scale 
dependence: 



^{L + L e )C VA ia 2 + 



A 2 



(10) 



with Ai, Ai constants (depending on the geometrical 
characteristics of the structure). Eq. (fT0|) shows that 
as the slabs length-scale becomes smaller the magnetic 
resonance frequency does not continuously increase, but 
beyond a length-scale it saturates to a constant value. 
The saturation of the magnetic resonance frequency is 
exactly what is observed in Fig. 2, and is exclusively due 
to the existence of kinetic inductance (note that the resis- 
tance, R, representing the ohmic losses, does not appear 
in the above formula (jT0|) ) . This inductance originates 
from the electrons inertia and represents the "difficulty" 
of electrons to follow high frequency motions, i.e. their 
difficulty to respond to high frequency fields. 

Note that without the consideration of L e the second 
term in the square root of Eq. (|10p would not exist, 
leading to a linear dependence u> ~ 1/a (as occurs in mi- 
crowaves and larger scales) making unable to explain the 
saturation behavior observed in Fig. 2. Note also that 
the kinetic inductance L e does not influence the response 
of the structure only at the magnetic resonance regime 
but at all high-frequency regimes, including permittivity 
resonances. 

Substituting Eqs. ((9]) in Eq. (|T0|) and defining nor- 
malized (dimensionless) geometrical parameters, i.e. V = 



1 1 a, t' = t/a, t' m = t m /a, one can obtain the dependence 
of ljlc from the geometrical characteristics of the struc- 
ture, as well as an expression for the saturation value, 



j lc- 



Bhl' 



e b V 2 



J LC- 



(11) 



A correction to the above formula (jTTJ) can be obtained 
if one takes into account also the potential energy of the 
electrons inside the metal slabs, through an equivalent 
capacitance^, C e = £owt m /l (the capacitance of the ca- 
pacitor formed inside the metal), added to the slabs ca- 
pacitance, C. In this case the saturation value for the 
magnetic resonance frequency is found as 



sat 



0J U 



ft' 



+ 1 



(12) 



showing that the absolute upper limit for the saturation 
frequency is not arbitrary high but it is restricted by the 
plasma frequency of the bulk metal. 

Since in the following we will consider systems with 
both V jt' and l'/t' m larger than unity, where the simpli- 
fied equation (fTTj) is still valid, we will omit the electron 
potential energy in the following discussion, keeping into 
account though that this energy /capacitance sets a finite 
upper limit for the saturation value of the magnetic res- 
onance frequency, which is the plasma frequency of the 
bulk metal. 



B. Magnetic permeability resonance 

As had been shown in Section II, by reducing the 
length-scale of the artificial magnetic structures the mag- 
netic permeability resonance becomes more and more 
weak, unable to lead to negative \i values beyond a length 
scale. This weakening is revealed in both the real part of 
/i (where smaller absolute values of the Max[Re(/^)] and 
Min[Re(^)] are observed) and the imaginary part (where 
smaller Max[Im(^)], at w = ujlc, is observed) - see Fig. 3. 

A detailed examination of Eq. ([7]) reveals that the 
strength of the resonance is determined from both fac- 
tors 7 and F'. On the other hand, the width of the neg- 
ative permeability regime, Aw, seems to be much more 
sensitive to the factor F' , and almost unaffected from 
7 (7 has only small influence on the lower limit of the 
negative permeability band, while F 1 strongly affects the 
upper limit of this band). Note that in the absence of 
losses, i.e. 7 = 0, the upper limit of the negative perme- 
ability band is ojlc/V^- — F'\ the lower limit is simply 
uj lc , i.e. 



1 



VI - F' 



1) 



(13) 
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Using Eqs. © and (0J in combination with with Eqs. 
((91), one can obtain the scaling dependence of both F' 
and 7, as 



F' = F- 



L 



L + L, 



cx BFa 2 



(B constant) and 



I? 



7 



L + L P 



1 

D ia 2 + D 2 



(14) 



(15) 



(Di, D2 constants). From Eq. (fT4)) one can derive two 
important conclusions: (a) The factor F', which mainly 
determines the frequency width of the negative perme- 
ability regime, is independent of the resistance, R, thus 
independent of any loss mechanisms, (b) With the con- 
sideration of the kinetic inductance the factor F' from 
scale independent (if L e is negligible) becomes scale de- 
pendent and tends to zero as the size of the structure 
becomes smaller and smaller (a — > 0). This means that 
even in the absence of ohmic losses, it would be impos- 
sible to get negative permeability values of band-width 
substantially larger than (a 2 / '\p)u> lc m arbitrarily small 
length scales. 

The scaling behavior of F' also implies that the relative 
bandwidth of the negative permeability regime, while it 
is maintained almost constant before the staring of the 
saturation (see Fig. 4), it becomes smaller and smaller 
deeper in the saturation regime, indicating that working 
before the saturation regime favors the widest bandwidth 
of each specific structure. 

The geometrical dependence of the factor F 1 for the 
slab-pair design can be easily obtained by substituting 
Eqs. m to Eq. 



F' = F- 



1 



1 



c 2 /(^ m i) 1 + A 2 /[(2tt) 2 ^ 



(16) 



Concerning the loss factor 7, from Eq. (fT5|) one can 
see that 7 increases as the length-scale decreases, justi- 
fying the higher losses in the smaller length scales. This 
increase though does not continue up to the smallest 
scales, but beyond a length-scale it approaches a satu- 
ration value. The geometrical dependence of 7 for the 
slab-pair design is obtained analogously with that of F', 



(17) 



1 + {Ul/C 2 )t m t 1 + {2TTft m t/\l 



Eq. p7]) shows that the saturation value of 7 is the colli- 
sion frequency of the bulk metal as considered in the free- 
electron (Drude) description of the metal^ 2 -. The satura- 
tion regime for 7 is approached simultaneously with the 
magnetic resonance frequency saturation, showing that 
deep in the saturation regime the weakening of the res- 
onance is not mainly the result of the ohmic losses but 



it is rather the result of the kinetic inductance (affecting 
through the factor F'). 

Finally, it is important to point out that both 7 and 
F' depend not only on the "quality" (i.e. plasma fre- 
quency and damping factor) of the metal used for the 
fabrication of metamaterials but also on the geometrical 
parameters of the structures (see Eqs. (Tj"6|) and l[i"?])). 
This reveals the possibility to modify these factors, and 
thus to enhance the performance of the high frequency 
metamaterials, by modifying the geometry. 



IV. OPTIMIZED DESIGNS 

From the analytical formulas and the discussion of the 
previous section it becomes clear that for achievement of 
optimized high frequency (e.g. optical) magnetic meta- 
materials one should require: (a) Highest possible sat- 
uration value for the magnetic resonance frequency; (b) 
larger possible parameter F', determining the strength 
of the magnetic resonance and the width of the negative 
permeability regime; (c) smallest possible loss factor 7. 

The general requirements for meeting the above con- 
ditions can be easily concluded based on Eqs. (p~0|) . (Ti"4j) 
and (|15p . They demand: Structures of small capacitance, 
C; structures of small kinetic inductance L e , compared to 
magnetic inductance, L; structures of low resistance, R. 
On the other hand, the role of the magnetic field induc- 
tance is more puzzling: while low inductance facilitates 
the achievement of high magnetic resonance frequency, 
it results to "lower-quality" resonance, i.e. weaker reso- 
nance and higher losses (compared to a higher inductance 
structure of the same length scale). Thus, the inductance 
optimization should be based on the specific requirements 
for the designed metamaterial. 

To translate the above optimization conditions to spe- 
cific geometrical and material conditions for the slab- 
pair- based systems we can use Eqs. (fTT|l . (fT6| and (jlTf . 
From these equations can be concluded that optimized 
slab-pair-bascd systems are favored from structures of 
thick metal (high t m ) and thick separation layer between 
the slabs of the pair (i.e. high t - not as high though as to 
cancel the interaction between the slabs) ; also structures 
of wide slabs (i.e. large width w - to maximize the struc- 
ture volume fraction F appearing in Eq. (flBj) V More- 
over, high quality optical metamaterials demand metals 
of the highest possible plasma frequency and the lowest 
possible collision frequency. 

The role of the slab-length is not one-way: while 
shorter slabs (compared to the corresponding unit cell 
side) facilitate high attainable magnetic resonance fre- 
quencies they lead to narrower negative permeability 
regime (due to reduced F in Eq. (fTl))) '). and vice versa. 

Finally, we should emphasize here that for enhanced 
metamaterial performance one should target operation 
below the saturation regime, as to ensure large negative 
permeability width and lower losses. 
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A. The fishnet design 

As was mentioned in the previous sections and can be 
easily concluded from Eq. (fTO]) . structures with reduced 
inductance lead to higher magnetic resonance frequen- 
cies. A structure based on slab-pairs which offers a con- 
siderable reduction in the inductance is the fishnet design 
- see Fig. 1(d). Studies of this design in microwaves^ 
have revealed that at the magnetic resonance loop cur- 
rents exist not only at the slab pair but also in the necks' 
part of the metallic element. This neck contribution can 
be taken into account in an effective LC circuit descrip- 
tion of the structure by considering an additional induc- 
tance, due to the necks, which is in parallel with the 
inductance of the slabs, resulting to a reduced total in- 
ductance, Lfishnet, with 



1 



1 



1 



-^fishnet L< 



la] >s 



knocks 



(18) 



(-^siabs and inocks represent the inductance of the slabs 
and the necks part respectively; note that the L s iabs here 
it is not exactly the same as that for the slab-pair-only 
system, due to the difference in the charge and current 
distribution between only-slab-pair systems and fishnet.) 
This reduced inductance leads to a magnetic resonance 
frequency 



1 



ui 



LC. fishnet 



L 



fishnet 



C 



J LC, slabs 



(1 



L 



slabs ■> 
-'necks 



(19) 



i.e. higher than that of only the slabs. 

By reducing the size of the structure down to submi- 
cron and nanometer scale, the inductance of both the 
slabs and the necks part gets a contribution from the ki- 
netic (electrons) inductance. This contribution though 
does not modify the general relation (|19j) . indicating the 
higher magnetic resonance frequency of the fishnet de- 
sign than that of the component widc-slab-pair system 
even in the saturation regime. 

This higher attainable magnetic resonance frequency 
makes the fishnet design the preferential one for the 
achievement of optical metamaterials, something that 
has been already revealed from many existing ex- 
perimental and theoretical works concerning optical 
met amat erial a 14 ' 15 . 

The relation (fTO)) . combined with Eqs. can easily 
lead to a relation for the geometrical dependence of the 
magnetic resonance frequency for the fishnet, i.e. 



J LC, fishnet 



J LC, slab 



s(l 



/ 



(20) 



and thus for the saturation value of this magnetic reso- 
nance frequency. In Eq. (f2U)l w n and w is the width of 
the neck and slab parts respectively (along H-direction) 
and <ze is the lattice constant along the E direction (see 
Fig. 1). Eq. (|2TH) suggests that for the achievement of 
high magnetic resonance frequency saturation values for 
fishnet, apart from the conditions for the optimization 



of the slab-pair components one should pursue also wide 
neck parts; moreover, the wider the neck parts the higher 
the saturation value of the fishnet magnetic resonance 
frequency is. 

Concerning the magnetic permeability expression for 
the fishnet, here the situation is more complicated com- 
pared to only-slabs systems, due to the more complicated 
current picture at the magnetic resonance^. Following 
the observations and conclusions of Ref<21, according to 
which the fishnet unit cell can be approximated with an 
RLC circuit with inductance of slabs and inductance of 
necks connected in parallel, the magnetic permeability 
for the fishnet can be calculated following the same steps 
as the ones presented in the previous section for the slab 
case, with modifications in the incident flux, the mag- 
netic moment and the total inductance per unit cell: 

Here the incident flux, </>, can be written as <p = 
/ioflE^oe^ 1 "*, and the modified magnetic moment as 



Islabslt— /necks(o.E~ l)t — -^fishnet [ 



It 



(oe - 0*1 



slabs 



-^nccks 

(21) 

^slabs-^slabs 



where / = I s i a bs + knocks and IL n 

shnet 
-^nccks -^nccks ■ 

The inductances L s iabs and L nec ks include both the 
magnetic field inductance and the electrons inductance 
for slabs and necks. 

From Eq. (|2T|) one can see already that the presence 
of the necks weakens the magnetic response of the struc- 
ture, since the neck contribution in the magnetic moment 
opposes that of the slabs. 

With the above considerations, the magnetic perme- 
ability for the fishnet design can be expressed as 



(l/V uc )Hoa E tuj 2 (lt/L s \ ahs - (a E - l)t/L necks ) ^ 

OJ 2 — 1/ifishnctC + twi?/L ns hnet 



fl((j) = /x [l- 

W- — 1/ ^fishnet O f IUJXI/ -Ivfishnet 

(22) 

Comparing the above equation with Eq. ([3]), one can 
observe that the factor multiplying the uj 2 in the nu- 
merator, which is the main factor determining the spec- 
tral width of the negative permeability regime, becomes 
smaller than that of the slab only case (F'). This can 
explain the reduced (compared to only slabs) spectral 
width of the fishnet negative /i regime observed in Fig. 
4. 



V. OBTAINING QUANTITATIVE RESULTS 

In the previous sections we presented a simple RLC 
circuit model for obtaining qualitative results and opti- 
mizing the slab-pair-based metamaterial structures. It is 
important to note though that the relations (3)-(4) do not 
only have a qualitative power, but they can be used also 
to obtain quantitative results, if plugged with more ac- 
curate expressions for the capacitances and inductances 
involved. 



10 



Here we will demonstrate this quantitative power of 
our relations in the case of the widc-slab-pair system. 
For wide-slab-pair systems (w = an) and examining 
the fields and currents at the magnetic resonance fre- 
quency, one can use for the total capacitance of the sys- 
tem the capacitance of two identical parallel plate ca- 
pacitors connected in-series, each of capacitance C = 
coeb(w/t)(l/2.5), and for the electrons inductance the 
same relation as in Eq. (9), multiplied by 2 (to take 
into account both slabs of the pair). 

The calculated magnetic resonance frequency obtained 
by using these relations is shown with the dotted-dashed 
line in Fig. 6, and as can be seen there it is in excel- 
lent agreement with the magnetic resonance frequency 
obtained through realistic calculations. 

Using the same equations for the case of narrow slabs, 
one can see that they fail to describe accurately the scal- 
ing of the magnetic resonance frequency. In fact these 
equations predict no dependence of w m on the width of 
the slabs, unlike the observed in Fig. 2 result, which 
shows lower magnetic resonance frequency for the nar- 
rower slabs. 

To explain this discrepancy, we should point out that 
as the slabs become narrower Eqs. (9) for the capaci- 
tance and the magnetic field inductance become less and 
less accurate, and in the limit of w < t they should 
be replaced by the corresponding expressions for par- 
allel wire systems, i.e. L = jjL^l ]n(t/w)/n, and C = 
e&eo7r(Z/2.5)/ln(£/'u;). The last relations for L and C 
lead to lower magnetic resonance frequency saturation 
values compared to the ones predicted by the correspond- 
ing Eqs. (9), in agreement with the numerical results 
presented in Fig. 2. 

VI. CONCLUSIONS 

In this paper we examine the magnetic response of res- 
onant magnetic structures based on the slab-pair design, 
as the structures are scaled down from mm to nm scale. 
This response is examined using detailed numerical sim- 
ulations for obtaining the frequency dependence of the 
magnetic permeability (including its resonant behavior) 
through reflection and tranmission data of realistic struc- 
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